Local defects and ferromagnetism in graphene layers. 



M. A. H. Vozmediano, 1 M. P. Lopez-Sancho, 2 T. Stauber, 2 and F. Guinea 2 ' 3 

1 Unidad Asociada CSIC-UC3M, Universidad Carlos III de Madrid, E-28911 Leganes, Madrid, Spain. 
Instituto de Ciencia de Materiales de Madrid, CSIC, Cantoblanco, E-28049 Madrid, Spain. 
3 Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, MA 02215, USA 

(Dated: February 2, 2008) 

We study the changes in the electronic structure induced by lattice defects in graphene planes. 
In many cases, lattice distortions give rise to localized states at the Fermi level. Electron-electron 
interactions lead to the existence of local moments. The RKKY interaction between these moments 
is always ferromagnetic, due to the semimetallic properties of graphene. 
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Introduction. A number of recent experiments sug- 
gest that pure graphite behaves as a highly correlated 
electron system[]J. In particular it shows a metal- 
insulator transition in magnetic fields and insulating be- 
havior in the direction perpendicular to the planes in 
different samples 0, 0. [J EL.0 111- Recent re- 
sults show ferromagnetic behavior |9j, enhanced by pro- 
ton bombardment [10|, what opens up a new way to 
the creation of organic magnets 11]. In this paper we 
study the formation of local moments near simple de- 
fects. It is shown that many tyoes of lattice distortions, 
like cracks and vacancies, can induce localized states at 
the Fermi level, leading to the existence of local mo- 
ments. The RKKY interaction between these moments is 
always ferromagnetic due to the semimetallic properties 
of graphite. Hence, the RKKY interaction does not lead 
to frustration and spin glass features. 

The model. The conduction band of graphite is well 
described by a tight binding model which includes the tt 
orbitals which are perpendicular to the graphite planes at 
each C atomE2- If the- interplane hopping is neglected, 
this model describes a semimetal, with zero density of 
states at the Fermi energy, and where the Fermi surface 
is reduced to two incquivalcnt K-points located at the 
corner of the hexagonal Brillouin Zone. The low-energy 
excitations with momenta in the vicinity of the Fermi 
points have a linear dispersion and can be described by 
a continuum model which reduces to the Dirac equation 
in two dimensions^! El El El 13 . The Hamiltonian 
density of the system is 

Ho = hv F J d 2 r$(r)(ia x d x + ia y d y )9(r) , (I) 

where the components of the two-dimensional spinor 
^(r) = (\&i(f), ^(f)) correspond to the Bloch states 
of the two sublattices in the honeycomb structure, in- 
dependent in the absence of interactions. In the clean 
two-dimensional system there is no room for low energy 
electronic instabilities, the short range interactions being 
irrelevant due to the vanishing density of states at the 
Fermi level. 

It is known that disorder significantly changes 
the states described by the two dimensional Dirac 
equationpjl El El, 

and, usually, the density of states 



at low energies is increased. Lattice defects, such as pen- 
tagons and heptagons, or dislocations, can be included in 
the continuum model by means of a non-abelian gauge 
fieldEl IH that reproduces the effects of the curva- 
ture of the lattice and the possible exchange of Fermi 
points. Within the same theoretical scheme it has also 
been shown that certain types of disorder randomly dis- 
tributed in the gra phene lattice enhances the effect of 
the interactions 2 lj and can stabilize new phases. In 
addition, a grap hene plane can show states localized at 
interfaces (23, [23, which, in the absence of other types of 
disorder lie at the Fermi energy. Structures with mixed 
sp 2 and sp 3 bonding can also lead to localized statespij. 

The tight binding model defined by the tt orbitals at 
the lattice sites can have edge states when the sites at the 
edge belong all to the same sublattice[22l I23I Elf (zigzag 
edge). These states lie at zero energy which for neutral 
graphene planes correspond to the Fermi energy. In the 
continuum model described earlier localized states are 
normalizable solutions (\E'i(f ), \E , 2 (r)) of the Dirac equa- 
tion: 

(d x + id y )*i(r) = 
(d x -id y )*2(r) = (2) 

These equations are satisfied if ^iX?) is an analytic func- 
tion of z = x + iy and ^(f) = 0, or if v E , r(f) = and 
^(f) is an analytic function of z = x — iy. Zigzag edge 
states can be obtained as follows. Consider a semi-infinite 
honeycomb lattice with an edge at y = and which oc- 
cupies the half plane x > 0. A possible solution which 
decays as x — * 00 is 

Vi(x,y) oc e~ kz = e lkv e- kx ,y 2 (r) = . 

These solutions satisfy the boundary conditions at y = 
if the last column of carbon atoms belong to the sub- 
lattice where the component ^1 is defined. Then, the 
next column belongs to the other sublattice, where the 
amplitude of the state is, by construction, zero. 

This kind of solutions can be generalized to describe 
other types of extended defects that will be produced in 
experiments where graphite samples are bombarded by 
protons. In a strongly disordered sample, large defects 



FIG. 1: Elongated crack in the honeycomb structure. The 
crack is such that the sites in the upper edge belong to one 
sublattice, while those at the lower edge belong to the other. 
Bottom: approximate cut in the complex plane which can be 
used to represent this crack at long distances. 



made up of many vacancies can exist. These defects will 
give rise to localized states, when the termination at the 
edges is locally similar to the surfaces discussed above. 
Note that, if the bonds at the edges are saturated by 
bonding to other elements, like hydrogen, the states at 
these sites are removed from the Fermi energy, but a sim- 
ilar boundary problem arises for the remaining tt orbitals. 
The only possible localized states can exist at zero en- 
ergy, where the density of extended states vanishes. The 
wave functions obtained from the Dirac equations will 
be normalizable and analytic functions of the variables 
z = x + iy or z = x — iy of the form 

= 0] 

obeying the boundary conditions imposed by the shape 
of the defect. 

Extended vacancies with approximate circular shape 
can support solutions of the type 

f(z) oc z~ n , n > 1 . 

By using conformal mapping techniques, solutions can be 
found with the boundary conditions appropriate to the 
shapes of different defects. 

A simple case is the elongated crack depicted in Fig.pp, 
which we assume to extend from i = -a to i = a, and 
to have a width comparable to the lattice constant along 
the y axis. The analytic function f(z) associated with 
localized states near a crack of this shape should satisfy 
Re/(z) = at the crack edges, because the boundaries of 
the crack include atoms from the two sublattices. Hence, 
the boundary of the crack leads to a branch cut in the 
complex function f(z). Labelling edge states by a quan- 
tum number n, we find that the function 'F can be written 
for these states as: 



* „ = <^ Re 



A 



z n \J z 2 — a 2 



,0 



FIG. 2: Energy of the five states with positive energy closest 
to e = for a cluster of size Lx L which contains at the center 
the defect shown in the inset. Periodic boundary conditions 
are used. The spectrum is symmetric around e = 0. 



A similar solution is obtained by exchanging the upper 
by the lower spinor component, and replacing z <-> z. 
Because of the discreteness of the lattice, the allowed 
values of n should be smaller than the number of lattice 
units spanned by the crack. 

We have checked numerically the existence of these lo- 
calized states by diagonalizing the tight binding hamilto- 
nian in finite lattices of different sizes. The dependence 
of some of the states close to the chemical potential (zero 
energy) on the cluster size is shown in Fig. [2]. The delo- 
calized states show a dependence £dci oc L , consistent 
with the properties of the Dirac equation from which they 
can be derived. The states closest to e = 0, show a depen- 
dence £i oc cx £~ 2 , which suggest a power law localization, 
in agreement with the previous analysis. 

These states are half filled in a neutral graphene plane. 
In the absence of electron interactions, this leads to a 
large degeneracy in the ground state. A finite local re- 
pulsion will tend to induce a ferromagnetic alignment of 
the electrons occupying these states, as in similar cases 
with degenerate bands 26] . 

To verify this behavior and to compare with previous 
results we have performed an unrestricted Hartree Fock 
calculation of the graphene cluster with different types 
of defects including large vacancies as the one shown 
in Fig. [p. Similar calculations done in graphene rib- 
bons have shown non-bonding molecular orbitals local- 
ized mainly along the zigzag edges (edge states) [23, |28| . 
These edges have partly flat bands which give rise to a 
sharp peak in the density of states at the Fermi energy. A 
finite on-site Coulomb interaction leads to an instability 
of these flat bands and to the existence of spin-polarized 
states. We have verified that the similar behavior can 
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FIG. 3: Density of states of a 26 x 26 cluster with a large 
defect with zig-zag boundary as the one shown in Fig.0 The 
on-site interaction term is U = 0. 
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FIG. 4: Same as Fig. [2] with an on-site interaction U = 0.5t. 



also be obtained near the boundary of a large defect as 
the one shown in Fig. [P . We have performed calculations 
in a honeycomb lattice of various sizes with a maximum 
of 26x26 unit cells, i. e. a cluster of 52 x 52 sites. In 
the absence of defects the density of states is zero at the 
Fermi level. Figure [3] shows the charge density distri- 
bution of a lattice with an extended vacancy shaped as 
shown in Fig.pQ, without on-site interaction. The Fermi 
level is located at zero. A peak in the Fermi level is 
formed due to the presence of the defect. The total mag- 
netization of the cluster is zero. When the on-site in- 
teraction is switched on, a spin polarization arises of a 
magnitude proportional to the size of the defect. Fig.0] 
shows the density profile of the polarized cluster obtained 
with U = 0M. 

Next we will analyze the influence of these magnetic 
moments on the conduction band. The hopping between 
the states involved in the formation of these moments and 
the delocalized states in the conduction band vanishes by 
definition, if the localized states lie at zero energy. Hence, 
a Kondo-like coupling mediated by the hopping will not 
be induced. The localized electrons act as a reservoir of 
localized moments, which interact with the valence elec- 



trons via the on site Hubbard repulsion, U. The wave- 
function of the localized electrons overlaps on many sites 
with the valence electrons. On these sites, the valence 
electrons will tend to be polarized with their spins paral- 
lel to the spins of the localized electrons. The change in 
the energy at a given site for the valence electrons near 
the defect is of order U, distributed over N sites. Then, 
the valence electrons will induce long range RKKY in- 
teractions between the localized moments, which can be 
estimated by adding the contributions from all sites. We 
find: 



Jrkky(v) ~ U 2 N 2 I riVMk) ~ U 2 N 2 



— Fi3 ( 3 ) 



Where a is the lattice constant, and the static suscep- 
tibility is given by: x(k) oc |k| 16], and a is the lattice 
constant. 

Due to the absence of a finite Fermi surface, the RKKY 
interaction in eq.© does not have oscillations. Hence, 
there are no competing ferro- and antiferromagnetic cou- 
plings, and the magnetic moments will tend to be fer- 
romagnetically aligned. The total polarization per unit 
area at low temperatures is proportional to Nc, where c 
is the concentration of defects, and N is proportional to 
the average size. 

We can make an estimate of the Curie temperature 
from the coupling between the local magnetic moments 
given in eq. . The entropy cost of aligning ferromag- 
netically moments is S ~ T per moment. The average 
distance between moments is |f| ~ c -1 / 2 . Hence, the 
free energy per moment in the ferromagnetic phase can 
be written as: 



U 2 N 2 a 4 c^ 2 . , 

F(m) = -a + c 2 T ] to 2 



(4) 



where c\ and c 2 are numerical constants of order unity. 
The value of the free energy will be negative (and below 
the value in the paramagnetic phase) at a Curie temper- 
ature given by: 



U 2 N 2 a 4 c 3 ' 2 _ U 2 Na 3 
~ F ~ ~W^ 3 ' 



(5) 



where W is the conduction electron bandwidth, W ~ 
Vf/cl, and I is the average distance between impurities. 

The Curie temperature depends on the concentration 
and size of the defects. Assuming, as an example, that 
N ~ 10 and / ~ 10 2 a, we obtain a saturation magnetiza- 
tion of 10 -3 Bohr magnetons per unit cell, and a Curie 
temparature T c ~ 10~ 4 U 2 /W. The value of U 2 /W can 
be estimated to be ~ leV. Then, these arguments give 
Tc ~ IK. This temperature is considerably lower than 
the experimentally observed ones. It is worth noting, 
however, that this analysis does not take into account 
the enhancement of the susceptibility of the conduction 
electrons, percolation effects due to the random distribu- 
tion of impurities, and the finite extension of the localized 
states induced by the defects. 
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Conclusions. We have shown that, under very gen- 
eral circumstances, lattice defects, vacancies and voids 
in the graphene structure give rise to localized states at 
the Fermi energy. The number of these states scales 
roughly with the perimeter of the defect. Repulsive 
electron-electron interactions lead to the polarization of 
these states, and to the formation of local moments. The 
RKKY interaction mediated by the valence electrons de- 
cays as r~ 3 , where r is the distance between defects, and 
shows no oscillations, due to the vanishing of the Fermi 
surface in a graphene layer. The interaction is ferromag- 
netic, and the system cannot show the frustration effects 
and spin glass features observed in other disordered sys- 



tems with local moments. On the other hand, the Curie 
temperature estimated assuming a random distribution 
of local moments is low, Tq ~ IK, for reasonable values 
of the defect concentration. It may happen that perco- 
lation effects, and the finite extension of the localized 
states which give rise to the local moments will increase 
the value of To 
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